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X ■ Abstract 

M'. 

We apply a recent theoretical analysis of hadronic observables in inclusive 
semileptonic heavy hadron decays to the phenomenology of B and D mesons. 
Correlated bounds on the nonperturbative parameters A and Ai are derived 
by considering data from B decays and, independently, data from D decays. 
The two sets of bounds are found to be consistent with each other. The data 
from B decays are then are used to extract a lower limit on the CKM matrix 
element \Vcb\- We address the issue of the convergence of the perturbative 
expansions used in the analysis, and compare our bounds on A and Ai to 
lattice and QCD sum rule results. Finally, we argue that a comparison of 
the analyses of D and Ds decays provides evidence for the applicability of 
parton-hadron duality in the semileptonic decay of charmed hadrons. 
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I. INTRODUCTION 



The heavy quark hmit of QCD is of enormous practical use, because with it one may 
describe a wide variety of heavy hadron decay rates and matrix elements in terms of a 
small number of parameters. These parameters reflect nonperturbative QCD effects and 
cannot be computed directly. Instead, they must either be modeled or, preferably, be ex- 
tracted from experimental data. One of the most important applications of the analysis 
of inclusive decays is the determination of the CKM matrix element |Vc6| from the process 
B Xcii', which is complementary to the extraction from the exclusive decay B — > D*tv. 
The computation involves an expansion in powers of l/mb, and to 0{l/ml) there appear 
three nonperturbative parameters: A (or equivalently, the quark mass rrih), Ai, and A2. 
While A2 may be extracted directly from the B-B* mass splitting, A and Ai are not directly 
measurable. Two approaches currently popular in the literature are to employ various QCD 
sum rules to estimate these parameters, and to use an analysis of inclusive semileptonic D 
decay to fix one linear combination of them. However, each of these methods has a severe 
disadvantage: the QCD sum rule results are not model-independent consequences of QCD, 
and the expansions in as{mc) and l/rric may or may not work well at the low scales relevant 
for D decay [|,g. 

In a recent analysis 0, we calculated the leading perturbative and nonperturbative 
contributions to moments of the hadronic energy and invariant mass spectra in semileptonic 
heavy hadron decays. These predictions are particularly interesting, because experimental 
information on invariant masses of the hadrons produced in these decays may be derived 
from the reported branching ratios to exclusive final states. Furthermore, they rest on the 
same theoretical basis as earlier analyses of other semileptonic quantities such as the decay 
rate and the lepton energy spectrum [^]. 

In Ref. we performed some preliminary phenomenology based on this theoretical anal- 
ysis, deriving correlated bounds on the nonperturbative parameters A and Ai. In this paper 
we will develop this phenomenology further, incorporating additional data and including in 
the discussion the semileptonic decays of charmed mesons. Our main conclusions are: 

1. The perturbation series appearing in the analysis are under better control than had 
previously been thought. While the two-loop corrections relating A and Ai to the 
semileptonic decay rate and to the first moment of the invariant mass spectrum in 
B —>■ XcCu are large, these corrections partially cancel in the relation between the 
semileptonic decay rate and the first moment of the invariant mass spectrum. The 
corresponding perturbation series relating the two physical quantities appears to be 
better behaved. Using the scale setting technique of Brodsky, Lepage and Mackenzie 
[^], we find a BLM scale /iblm = 0.26 for the relation between the first moment 
and the semileptonic width. This extends our previous result to the case of finite 
charm quark mass. 

2. When combined with the measured semileptonic width of the B, the moments of the 
invariant mass spectrum and the measured branching fraction to excited states yield 
the constraint 




(1.1) 
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While this is consistent with previous determinations of \Vcb\ from inclusive B de- 
cays , this result differs from previous extractions in that it does not depend on any 
assumptions about the size of A, nor on QCD sum rule estimates of the quark masses. 

3. The values of A and Ai extracted from the semileptonic decay width and first moment 
for D and decays are consistent with those obtained from B decays. The B results 
are also consistent with recent lattice extractions of the MS mass fni){mb) 0,^]. The 
combined results from B and D decays are inconsistent with the large negative value 
of Ai extracted from certain QCD sum rules 0. 

4. The theoretical prediction for the difference of the first moments of the invariant 
mass spectrum in D and Dg decays is well behaved and provides a test of parton- 
hadron duality. The comparison with experimental data is quite successful, providing 
additional evidence for the applicability of duality to the decays of charmed hadrons. 

Since our conclusions on D decays disagree significantly with those presented in Ref . , 
it is worth commenting on the discrepancy. The authors of Ref. used the extraction of 
rrib from Ref. [IlOi along with the QCD sum rules extraction of Ai = —0.6 ±0.1 GeV^ M 



to conclude that 1.25 < mP°^^ < 1.40. This results in a semileptonic decay width for the 
D meson which is at least a factor of two smaller than observed. However, it is difficult to 
relate this extraction of the pole mass to physical quantities. The radiative corrections in the 
relation between mj?"'*^ and the semileptonic charm width are so large that the perturbation 
series appears uncontrolled; whether or not this is the case for the relation between the 
moments of a{e'^e~ hh) (from which m^"^^ , and hence m^"'*', are extracted) and the 
charm quark semileptonic width requires a higher order calculation. Given this uncertainty, 
we prefer to treat Ai and A as free parameters, to be fixed by relations between the decay 
widths, moments and MS masses, in which case we find that all the data on charm and 
bottom are consistent with the smaller value Ai ^ —0.1 GeV^. This is also consistent with 



the observations of Ref. |rT[, where it is argued that the correct QCD sum rule should give 
a substantially smaller value of Ai than that found in Ref. 0. 

Finally, we note that we will consider values of Ai which violate the constraint Ai < 



-3A2 ^ -0.36 GeV^ which was proposed in Ref. [|12[. The proof of this bound is criticized 



in Ref. [|T^], where it is demonstrated that the inclusion of radiative corrections precludes 
any rigorous constraint on Ai. Hence we do not apply this proposed limit to our analysis. 



II. CONSTRAINTS FROM B DECAYS 



A. Theoretical expressions 



We begin by discussing the constraints which may be obtained from inclusive semileptonic 
B decays. The theoretical treatment of these decays involves a double expansion in powers 
of asifTib) and l/m^, employing an Operator Product Expansion (OPE) and heavy quark 
symmetry. From Ref. we have the expressions for the first two moments of the hadronic 
invariant mass spectrum for the process B Xciv, 
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where = 0, and we have defined the spin- averaged D and B meson masses, 
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In deriving the expressions ( |2.1|) , we have ehminated the ratio of pole masses mc/m^ by 
instead writing the heavy quark expansion in terms of mD/mB, 
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Performing a similar substitution in the expression for the semileptonic decay rate, we find 
(for m^ = 0) 
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The advantage of writing rs.i.(-B) and the moments ( |2.1| ) in this way is that there is now no 
hidden dependence on the heavy quark masses; the coefficients arising at each order in the 
OPE are determined by measurable quantities. 

The moments of the invariant mass spectrum depend only on the nonperturbative pa- 
rameters A, Ai and A2, and on the strong coupling constant aJmb) at leading order. Since 
^2{'n^b) = 0.12 GeV^ is known from the B-B* mass splitting,lil and as{mb) is measured in 
other processes, these moments provide direct information on the unknown hadronic matrix 
elements A and Ai. This information may then be inserted into the expression for rs.i.(-B) 
to determine \Vcb\ from the measured decay rate. 



B. BLM scale setting for finite nic 

We begin by addressing the question of whether the perturbative corrections to the 
relation between the semileptonic decay width and the moments of the hadronic invariant 



^In this paper, we will neglect the small running of )^2{p) between rub and rric. 
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mass spectrum are in fact well behaved. Earlier analyses have indicated that the two-loop 
corrections to Fg i (5) are uncomfortably large. In these analyses, one computes that part 



of the two-loop correction which is proportional to the first coefficient Pq = 11 



in the 



QCD beta function, and from this derives a BLM scale 
result 



for the process. One finds the 
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which, since /5oas(m6)/7r ~ 0.6, leads to a perturbation series which is quite poorly behaved. 
Following the BLM prescription of absorbing this term into the 0(as) correction by a change 
of scale, one finds a low BLM scale, /xblm = 0.16mb ^ SOOMeV.cl 

In Ref. 1^ we discussed a similar situation in the analysis of the decay b — > uiu. There 
we considered two perturbation series, neither of which is particularly well behaved: 
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The BLM scale for r(_B XJu) is /xblm = 0.08 rrib, while for {sh) it is /xblm = 0.03 mft. 
However, both of the expressions ( |2.6|) depend on the nonperturbative parameter A, which 
is defined only up to certain arbitrary conventions [111]. If the poor convergence of the 



perturbation series can be absorbed into A, then the large higher-order terms will be of no 
consequence, since ultimately A is eliminated from relations between physical observables. 

In Ref. Q we investigated whether this might be so by eliminating A from the equations 
(^.61), solving for V{B — > XJu) in terms of {sh)- Doing so, we found 
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leading to a much higher BLM scale, /Ublm = 0.38 m^. The apparent convergence of the 
perturbation series improves considerably under such a reorganization. 

For finite charm quark mass we may perform a similar analysis. We use standard tech- 
niques [15| to extract the two-loop term of the form /5o(as/7r)^ which contributes to the first 
moment {sh — tt^'d)- The calculation is straightforward but tedious, with the final integrals 
performed numerically. We find 



^This scale arises from taking mc/mi, = 0.37 (see Eq. (|27 
of Ref. IQ, where mc/rrib was taken to be 0.3. 



and differs slightly from the result 
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which again leads to a perturbation series which appears to be badly behaved, with a very 
low BLM scale, /isLM = 0.02 m^. However, if instead we use the expression ( pi^ ) to eliminate 
A from the semileptonic width (p.5|), we obtain 



XJu) 



Gl\V,t\^m% 
1927r3 



0.369 



1.17 



-0.74 



TT 



TT 



m. 



(2.9) 



The two-loop correction in Eq. ( |2.9| ) has been reduced by a factor of almost two compared 
with that in Eq. (|2.5|) . The one- loop correction is reduced as well, so the change in the 
BLM scale is less dramatic; the new BLM scale is ^blm = 0.28 m;,. This reorganization 
of the perturbation series gives us hope that the expansion a;s(m^)/7r is now under control, 
although, of course, the full 0{a1) correction remains an important source of uncertainty. 

For the second moment of the invariant mass spectrum there is, as in the massless case, 
no such cancelation. We find, for ^^(mb) = 0.22, 
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Since the 0{K) term comes with an explicit factor of a^, substituting a physical quantity for 
A will not introduce a term of (9(a^/?o) to cancel the large two-loop correction in Eq. ( p.lOj) . 
Therefore, we expect that constraints from the second moment will be more sensitive to 
higher order corrections than those from the first moment and hence less reliable. Fortu- 
nately, the most useful constraints in the A — Ai plane will come from the first moment of 
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C. Bounds on A and Ai 

Using the theoretical expressions ( |2.1| ) and experimental data, we now derive constraints 
on the nonperturbative parameters A and Ai. These quantities dependent on the scheme by 
which perturbation theory is defined; the bounds which we will derive are for A and Ai at 
one loop in QCD in the MS scheme, with the renormalization scale /i = m^. We make no 
claim that this is the "natural" definition of these quantities, and in any case the scheme 
dependence drops out of relations between physical observables. However, although they are 
unphysical, it is convenient to retain these parameters in intermediate stages of calculations, 
and in order to compare the values of A and Ai obtained from different observables we must 
specify some convention for their definition. 

In Ref. 0, we used the known branching ratio of B mesons to excited charmed mesons 
to estimate experimental lower bounds on {{sh — fn'jj)^). This estimate was based on the 
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OPAL measurement [0 of 34 ± 7% for the fraction of semileptonic decays to the states 
Di and However, while the sum of the two branching fractions is consistent with the 
recent CLEO 90% c.l. upper hmit |T^ of 30%, there appears to be a discrepancy with the 



branching fractions to the individual Di and D2 final states. In Ref. we took the average 
invariant mass of the produced Di(2420) and 1^2(2460) states to be 2450 GeV. Here we will 
assume that only the lower mass Di is produced, giving a more conservative lower limit on 
{{sh — ^i))")- We will also take the la OPAL lower limit on the fraction of semileptonic 
decays, 27%, so as to be consistent with the CLEO result. Doing so, and using the results of 
Ref. 1^ for the contribution to the moments from the D and D*, we find the experimental 
lower limits 

{sh - m^)min = 0.49 GeV^ , 
{{sH -mlf)min = l.lGeV\ (2.11) 

Note that in obtaining these limits we have assumed that no other excited states are pro- 
duced. It is more realistic to assume that there will also be production of the p wave doublet 
Dq and Di, which will raise the average invariant mass of the final hadronic state. However, 
since there is no experimental information on these states, we are conservative and do not 
include them in our estimates of {{sh — fn1)"')jain- 

Another observable which depends on A and Ai is the ratio of partial widths 

rjB^X^ (2,12) 

The theoretical expression for Rr also depends on the ratio of masses rrir/rnh, both at tree 
level and in the nonperturbative |r8| , pJ[] and perturbative p0|j21[] corrections. As before, 
irir/mi, may be re-expanded in terms of the observable rrir/mB, yielding the result 
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At present, there are only data on the average h hadron semitauonic branching fraction, ob- 
tained at LEP, where the identity of the bottom hadron is not determined. The experimental 
result is 

Sr(6^Xcrz/) = 2.75 ±0.48%. (2.14) 

This differs from Br{B —>■ XcTv) by contamination from the Eg and A^. However, this 
difference is small compared with the experimental uncertainty. The ratio of the theoretical 
expressions for R^^ in the B and A;, sectors is 

1 + (2.15) 



R;r' \ m 



B 



Since only about 10% of b hadrons at LEP are A^'s, the effect on R^ should be much less 
than 1%. Hence we use the measurement ( p.l4|) , along with Br (5 X^ev) = 10.7 ± 0.5 



24] , to obtain 

i?^ = 0.26 ±0.05. (2.16) 
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FIG. 1. The limits on A and Ai from data on semileptonic B decays. The sohd curves are lower 
limits on A from the first two moments of the hadronic invariant mass spectrum, while the dashed 
curve is a la upper limit from the ratio Rr- The solid curve on the left corresponds to the bound 
from {sh — Tn'jj), the one on the right to the bound from {{sh — "^d)^)- 



The comparison of the theoretical predictions (|2.1|) and ( p.l3|) with experiment leads to 
limits on A and Ai. The experimental central value for R.^ yields a curve which is entirely 
inconsistent with the other data, giving a negative value for A. Therefore in Fig. || we show 
the curve corresponding to the la lower limit on R^, along with the constraints from the 
moments of the invariant mass spectrum, where we have taken as{mb) = 0.22. Since the 
experimental error on R^ is relatively large, the 2a constraint is uninteresting, allowing all 
values of A and Ai in the displayed region of Fig. |I[ Therefore, at present we can only 
conclude that Rr favours a negative value of Ai. However, if the experimental uncertainty 
in R-r is reduced in the future, it may become an important quantity for constraining A and 
Ai. For now, the most interesting constraints in the A — Ai plane come from {sh — ^"d)- 

We note that a very similar discussion of the limits on A and Ai which may be obtained 
from R-r has been given by Ligeti and Nir ||23|. Our analysis is organized somewhat differ- 



ently from theirs in its treatment of experimental masses and errors, leading to results of a 
superficially different form, but the physics, and the uncertainties, are largely the same. 



D. Constraints on \Vcb\ 



Leaving aside the weak constraints from R^, we now take the information on A and Ai 
obtained from the analysis of the moments of sh and apply it to the extraction of \Vcb\ 
from the semileptonic width. We use the theoretical expression ( |2.4| ), the experimental 
semileptonic branching ratio of 10.7% p2|, the central value tb = 1.60 ps for the B lifetime, 
and the strong coupling constant as{mb) = 0.22. The experimental lower limit on (sh—^'d) 
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gives a restriction on the one-loop value of A 
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Including the 0{a'j.Po) term and assuming that this dominates the two-loop result, we obtain 
for the two-loop value of A the constraint 
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Incorporating the latter bound into our expression for the semileptonic width and solving 
for iV^bl, we find 
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In the second hne above we display the tree-level, one-loop and (partial) two-loop contribu- 
tions to the bound. As we showed earlier (see Eq. (^^) ), the perturbation series appears to 
be well behaved. We also note that this lower limit on \Vcb\ is relatively insensitive to the 
experimental error on {sh — ttt^d)- semileptonic production of Di and D2 mesons is 

reduced to its 2a OPAL lower limit of 20%, then the bound on A is weakened by 100 MeV. 
However, we still obtain > [0-038 - 2.8 x 10-^(Ai/0.1 GeV^)](rB/1.60ps)-i/2_ 



E. Constraints on mfe(m;,) 

Our bounds on A and Ai may be translated into bounds on the MS quark mass mb(mfe). 
However, these results should be treated with some caution because the two-loop corrections 
between mb(mft) and the first two moments oi sh — ^"d quite large, indicating a poorly- 
behaved perturbation series. 

We define the dimensionless parameter 

niB - mbirrib) A 4 a,(mb) , , as(^fe) V , .r, r,n\ 

xb = = h 77 h 1.56/5o +..., (2.20) 

m-B f^B O IT \ TT J 

where we keep only the large O (a^/3o) contribution to the two-loop term. We then find 
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i^E^^ = ( 1 - 0.52^ + l.Uxs] - 0.26^ - 0.26/5o ( 
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+ 1.13^ + 0.03^ + ... (2.21) 
for the first moment, and 

2 = —0.11 h O.Oooa;^ + 0.0053 

+ 0.0078/?o ( ^^^] ' - 0.14^ + . . . (2.22) 
\ TT J mj^ 

for the second. The O{a1(3o) corrections are clearly substantial. The corresponding BLM 
scales in Eqs. ( p.21| ) and ( p. 221 ) are /iBLM = 0.14 and 0.05 m;,, respectively, corresponding 
to O{all3o) terms which are roughly 60% and 80% of the one-loop term. Therefore there 
are likely to be much larger uncalculated radiative corrections in the relations between the 
first moment and rnb{mfy) then between the first moment and the semileptonic B width. 

With these caveats in mind, we combine these results with the experimental limits (|2.11|) 
to yield the constraints on rnh[mh) and Ai shown in Fig. |^. To illustrate some of the remaining 
dependence on the renormalization scale fj,, the constraints are plotted for both fj, = rub and 
fi = mf,/2. We also display, with the quoted uncertainties, the results of two recent lattice 
extractions of mi,{mb), 

7rif,(mb) = 4.17 ±0.06 GeV § , (2.23) 

and 



mb{mb) = 4.0 ± 0. 1 GeV || . (2.24) 

We see that our bounds are consistent with these lattice results. However, it is much more 
difficult to compare these limits to the extraction of the b quark pole mass from high moments 
of a(e+e- ^bb) fT^, 

= 4.827 ±0.007 GeV. (2.25) 

This is simply because the ambiguity in the renormalization scale introduces a large uncer- 
tainty, of order /3o(a<j/7r)^mfe ~ 200 MeV, in the relation of the "one-loop" pole mass ml"^'^ 
to mb{mb)- For example, if we use as{mb) in the one-loop relation 

mr'^^m,(m,)(l + ^^) , (2.26) 

then Eq. (|2.25| ) yields the value mb{mb) = 4.42 GeV, somewhat higher than the lattice 
results and barely consistent with our analysis of the moments. Furthermore, such a value 
certainly would be inconsistent with the combination of the moments analysis and the QCD 
sum rules extraction of Ai = —0.6 GeV^. On the other hand, in Ref. |jlO| it was argued that 



the natural scale for matching onto the non-relativistic theory is /i ~ 0.63 m^. Taking this 



10 



3.8 



-0.6 -0.4 -0.2 0.2 

Xi (GeV2) 

FIG. 2. The two-loop constraints on mf,(mfe) and Ai from the first two moments of sh — ^'d, 
evaluated at // = mj, (solid lines) and /i = Tnb/2 (dashed lines). The hatched region is excluded. 
Note that the constraints become more stringent as fj, is lowered. The two solid bands are the 
lattice results from Refs. [^,^ with the quoted uncertainties. The two-loop calculation includes 
only terms of order a^Po. 

lower scale, we find rnb{mb) = 4.34 GeV, still somewhat higher than the lattice calculations 
and still inconsistent with a large negative Ai. Of course, using an even lower scale in the 
one loop relation between m^°'^ and mb{mb) would lower the extracted value of mb{mb) even 
further. While this scale ambiguity is formally of higher order in than our calculation, 
we see that numerically it is quite significant. Without a higher loop calculation of m^"^'^ 
from QCD sum rules, it is impossible to determine whether or not this extraction of mf"^^ 
is consistent with the other constraints. 



III. D MESON DECAYS 

In Ref. it was argued that the semileptonic decays of charmed mesons are not well 
described in the heavy quark expansion, since the value of rric which is required to fit the 
observed semileptonic D decay rate lies significantly above 1.4 GeV. This is the upper limit 
suggested by combining the value of of rrib extracted from the T spectrum |l^ and the large 
negative value of Ai found from certain QCD sum rules 0. However, we believe that this 
argument should be reconsidered in light of the uncertainty inherent in relating the pole 
mass derived in Ref. [jlO[ with physical quantities. Indeed, we will find that the values of 
Ai and A implied by the semileptonic decay rate and the first moments of the D and Dg 
invariant mass spectra are in reasonable agreement with the limits from the corresponding 
observables in the bottom sector. 

The theoretical analysis in the D sector proceeds as before. Since is of order Aqcd, 
for consistency we keep only terms of order ml /ml in the theoretical expressions for (s^) 



and r{D Xgev). As we are neglecting terms of order OgAi and asA2, we also omit terms 
of order a^mg but keep logarithmically enhanced terms of order a^mg ln(m^/m^). Thus we 
find for the Cabibbo allowed semileptonic width, 
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A2 + Ai - A; 



(3.2) 



Note that the large infrared logarithms of the pole mass may be absorbed naturally into 
the MS mass renormalized at mc, fns{mc). As one would expect, the individual terms in the 
perturbative expansion, which arises from an operator product expansion performed at the 



scale ^ 



mr 



remain insensitive to physics below mc 



In addition to the usual nonperturbative parameters, the theoretical expressions for the 
decay rate and moments also depend on the strange quark mass ms{mc). We use the range 



100 MeV < m,(l GeV) < 300 MeV , 



(3.3) 



given by the Particle Data Group |25| (which changes only slightly when evolved from 
fi = 1 GeV to fi = mc). Finally, note that we have not re-expanded the m^ term appearing 
in the semileptonic width (|3.1| ) in terms of the meson mass m/). This is because here we 
have no analogue of the expansion ( p.3|) , since the strange quark is not heavy. As a result, 
the parameter A does not appear explicitly in Eq. ( p.l|) . Instead, we will solve Eq. ( p.l|) 
directly for mc, linearize in Ai, and then use the heavy quark expansion ( p.2|) to relate mc 
to A and Ai. 



The inclusive semielectronic D branching fraction recently has been measured to be |24 



Br{D^ Xe+u) = [6.64 ± 0.18{stat.) ± 0.29{syst. 



(3.4) 



The same CLEO analysis indicates that the Cabibbo-allowed inclusive rate is saturated, 
within errors, by the exclusive modes {K~ , K*~) e^u, with stringent upper limits 

reported on the channels —>■ (i^f (1270), {1A30)) e^v. The data on these decays yield 
for the first moment of the D — > Xgd'^i' invariant mass spectrum the value 



[sh)d = (0.49 ± 0.03) GeV^ 



(3.5) 



where we have added errors in quadrature and neglected the widths of the excited K mesons. 

The result of this analysis is the set of constraints displayed in Fig. |^, in which we 
also show the limits obtained earlier from B decays. Both theoretical and experimental 
uncertainties are included in the displayed bands. In fact, the dominant uncertainties are 
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FIG. 3. The restrictions on A and Ai from the semileptonic D decay rate (darker shaded region) 
and from {sh) in D decay (lighter shaded region). We also show the lower bounds on A from the 
moments of the B decay spectrum (cross-hatched black lines). 

theoretical, and have two distinct sources. First, there is the uncertainty in the strange 
quark mass (|3.3|) . Second, there is the effect of uncomputed terms in the mass expansion of 
order which will be more substantial than in bottom decays. The theoretical analysis 

at order is quite complex and involves a number of new nonperturbative parameters, 

so we do not attempt to include these terms systematically.! Instead, we obtain a minimal 
estimate of the size of the uncertainty arising from these effects by extracting the bounds 
from F(D — ^ Xgi'^u) in two ways: on the one hand, by solving directly for the width in terms 
of the charm quark pole mass, and on the other, by proceeding through the intermediate step 
of calculating the "decay mass" ^ . These two procedures, which are formally the same 
only up to order yield bounds on A which differ by approximately 70MeV. It would 

be hard to argue convincingly that effects were intrinsically smaller than this. For the 
analysis of {sh)d, we employ the simpler procedure of including a term n(0.5 GeV /mi))^ in 
the theoretical expression (|3.2|) and varying n between —1 and 1. 

Note that the constraints from the charm sector are compatible with those derived from 
bottom decays, although they do not appear to agree with the lattice extractions ofm;, 
shown in Fig. |^. However, Fig. |^ must be interpreted with caution, because of the possible 
presence of large higher loop effects. As discussed in Section IIC, each curve in Fig. ^ may 
be interpreted as giving a constraint on A derived from a particular process. A reliable 
comparison between the constraints derived from two different physical quantities requires 
that when one is expressed in terms of the other, the perturbative expansion in be well 
behaved. For example, we demonstrated earlier that although there are large two-loop 



^In Ref. Q the l/m|? corrections to the semileptonic decay rate were estimated using the factor- 
ization hypothesis and other arguments, and found to be small. 
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corrections to the theoretical expressions for Fs.i. and (sh) for B or D decays alone, these 
corrections partially cancel out when A is eliminated and one of these physical quantities is 
expressed in terms of the other. Unfortunately, this requirement may no longer be satisfied 
when expressions from the bottom and charm sectors are compared with each other. For 
example, writing {sh — 7n%) for B decays in terms of (sh) for D decays, we find 



{sh - ml,) = m 



B 



0.028 ^ ^ + 0.029/3o +0.12-^^^ 



TT \ VT / m 



(3.6) 



resulting in a pertubative expansion which does not appear to converge well. Thus while 
the rough consistency of the constraints from B and D decays is encouraging, it may not 
be particularly significant. Still, it is amusing to note that if we were to take the combined 
constraints as legitimate, then we would conclude that Ai is small and negative, of order 
—0.1 GeV^. Hence it would make a negligible contribution to most observables, including 

We do not compare the results from the charm sector with fnc{mc), which may be ex- 
tracted from the lattice measurement of fniy{mh) using the heavy quark mass relations. This 
is because the radiative corrections between these quantities are so large that perturbation 
theory appears uncontrolled, making it difficult to conclude whether or not the regions in 
the A — Ai plane indicated by the different observables are consistent with each other. 

A better test of duality in charm decays comes from comparing D and Dg decays. There 
is only an upper bound on the semileptonic branching ratio of the Dg 

Br{D, Xe+v) < 20% , (3.7) 

which is not strong enough to provide an interesting constraint from the semileptonic decay 
rate. The observed exclusive semileptonic modes are Ds [t], i]', 0) i'^u, for which CLEO 



has recently reported relative branching ratios |]26|. Adding the reported systematic and 
statistical errors in quadrature, we estimate a value for the first moment of the invariant 
mass spectrum, 

{sh)d^ = (0.68 ± 0.03) GeV^ . (3.8) 

We have assumed that decays to rj, rj' and (p saturate the inclusive semileptonic rate, an 
approximation for which we do not assign an error. We merely note that since semileptonic 
decays are known to be saturated by relatively few exclusive modes, the error due to 
this approximation for Dg might well be small. In light of this uncertainty, the upper limit 
on {sh)ds is perhaps somewhat less firm than the lower limit. 



The theoretical prediction for {sh)ds is given by the expression (|3l^) for {sh)d, with the 
replacements of m^i by mj)^ and A and Aj by their strange counterparts A^ and Aj^. They 
are related to A and A, by 

A. - A = ~ - ~ ^ 95 MeV , 

rriB - rriD 

Xis - Ai = [(ruB^ - fuB) - {mn^ - m^)] ^ -0.02 GeV^ , (3.9) 

ms - niD 

1 



A2s — A2 — - 



-0.01 GeV% 
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up to corrections to the mass expansions of relative order 1/mf ^. 

The theoretical expression for A{sh) = {sh)ds ~ {sh)d is particularly well behaved, 
since it vanishes up to SU{3) violating effects. Expanding in powers of l/rriD, we find 



A{sh) =mi)(A, -A) — + 




3 
10 



A) 



{\2s - A2) . 



(3.10) 



The large radiative corrections to (sh) almost precisely cancel in the difference. 

The leading corrections to this result likely come from the SU{3) violating Cabibbo- 
allowed annihilation channel in Dg decays. Although this channel typically has small sh, 
it contributes to (sh) Ds primarily through its effect on the total semileptonic rate. For 
small lepton mass, helicity conservation suppresses the purely leptonic decay Dg iv in 
favor of processes in which at least one gluon is emitted. Using vacuum saturation to guess 
the relevant hadronic matrix elements, one obtains a naive estimate of Airasirric) f^/m'j-) for 
the fractional correction to the Dg semileptonic decay widthi. With /^j ^ 200 MeV and 
as{mc) ~ 0.32, this corresponds to an increase of order 5% in the semileptonic width of the 
Ds, or a. decrease of {sh) by roughly the same percentage. Since the error in this estimate due 
to the use of vacuum saturation is unknown, we will take the effects of weak annihilation 
into account by decreasing the predicting for {sh)ds by 5%, and assigning an additional 
error of ±5%, or 100% of the naive estimate of the l/m\) correction, to the prediction of 
the differences in the moments. 

With the known values (|3^) of A^ — A and Aj^ — Aj, the expression (|3.10|) reduces to 
A{sh) = [0.136 + 0.121(A/1 GeV)] GeV^. For an analysis entirely within the charm system, 
it is appropriate to estimate A from the overlap of the shaded bands in Fig. 0, from which 
we obtain 0.25 GeV < A < 0.45 GeV. Including the estimated contribution and uncertainty 
from the annihilation channel in Dg decay, this range leads to A(sh) = (0.14 ± 0.04) GeV , 
which we may combine with the measured value of {sh)d to obtain the theoretical prediction 



This value agrees quite well with the experimental result (0.68 ± 0.03) GeV^. 

Indeed, we view the successful prediction of {sh)ds as evidence for the applicability 
of parton-hadron duality to inclusive semileptonic charm decays. The prediction of this 
moment is on a firmer theoretical footing than other quantities in the charm system, since 
the large radiative and power corrections up to 0{l/ml) cancel out of the difference of the 
moments. This feature distinguishes this test of duality from those in which a comparison is 
made to the bottom system, while its relative insensitivity to the value of A makes it more 
stringent than the comparison of the first moment to the total D decay width. At the same 
time, duality is satisfied in a nontrivial way, as the exclusive final states in D and Dg decays 
are entirely different. We are encouraged by this success of the theoretical analysis. 



^The authors of Ref. |Q] have calculated a related quantity of the same order, namely the contri- 
bution of "penguin- type" annihilation processes to the semileptonic decay of the D. They obtain 
the naive estimate above, multipled by a coefficient |(log(mc/^) + 1/3) ~ 1.3 for ~ 500 MeV. 




(3.11) 
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IV. SUMMARY 



In this paper we have explored the constraints on the nonperturbative parameters A and 
Ai which are obtained from semileptonic B and D decays. We have found that independent 
analyses of the bottom and charm systems yield limits which are consistent with one another. 
Taken together, they imply values of the order A ~ 450 MeV (at one loop) and Ai ~ 
-0.1 GeV^ Whether or not one chooses to trust the numerical results of the charm analysis, 
we see no evidence that parton-hadron duality fails in these decays. On the contrary, our 
discussion of the difference of the first moments in D and Dg decays leads us to quite the 
opposite conclusion: in at least one nontrivial case, duality works well for charm. 

A primary motivation for investigating inclusive decays is to extract the CKM matrix 

element |Vcfe| with high precision. Our analysis yields the lower limit \Vcb\ > [0.040 — 2.9 x 
10-^(Ai/0.1GeV2)](rs/1.60ps)-^/2 

using the current measurement of the branching fraction 
to excited D meson states in B decays. This is consistent with the value of \Vcb\ obtained 
from exclusive B decays 0. We have bolstered our theoretical analysis with a partial 
treatment of two-loop corrections to this bound, performing a BLM scale setting analysis 
which indicates that the relevant perturbation series is reasonably well behaved. 
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